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ABSTRACT 

We describe a population of small, high velocity, atomic hydrogen clouds, 
loops, and filaments found above and below the disk near the Galactic Center. 
The objects have a mean radius of 15 pc, velocity widths of ~ 14 km s _1 and 
are observed at \z\ heights up to 700 pc. The velocity distribution of the clouds 
shows no signature of Galactic rotation. We propose a scenario where the clouds 
are associated with an outflow from a central star-forming region at the Galactic 
Center. We discuss the clouds as entrained material traveling at ~ 200 km s _1 
in a Galactic wind. 

Subject headings: ISM: clouds — Galaxy: center 
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1. Introduction 

Many galaxies with active star formation in their central regions have large-scale Galac- 
tic winds. These are thought to be responsible for circulating enriched gas to the halo of 
their host galaxy and even to the intergalactic medium. The galaxy M82 has not only hot, 
ionized gas in its outflow, but also cool Ho; emitting filaments ( Bland fc Tully|1988 Shopbell 
fc Bland-Hawthorn||1998 ), and cold molecular gas (e.g. Nakai et aL||1987 Taylor et al.|2001 ). 
Simulations of starburst winds have shown that dense, molecular gas may be entrained in 
outflowing winds producing cold clouds and filaments at high latitudes around the center of 



a galaxy (e.g. Cooper et al. 2008 Melioli et al. 2013). 
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There is increasing evidence that the Milky Way may also have a large-scale nuclear 



wind. Bland-Hawthorn & Cohen (2003) linked the Galactic Center lobe and North Polar 

2° of the Galaxy. In 
~ 60°, centered on the 



Spur to a bipolar wind traced by infrared emission in the inner ^ 
gamma-ray emission there are two very large lobes extending to \b\ 
Galactic Center ( Su et al.|2010 ). Recently Carretti et al. (2013) identified counterparts to the 
Fermi bubbles in radio polarized emission and suggested that the lobes are formed through 
the effects of massive star formation, which drives a Galactic wind at ~ 1100 km s _1 . There 
are few direct tracers of the velocity of the putative nuclear wind of the Milky Way, and little 
evidence of associated cold gas. In fact, the inner 3 kpc of the Milky Way shows a distinct 
lack of diffuse H I at \z\ < 400 pc, which Lockman (1984) attributed to a Galactic wind 



removing H I except in the narrow disk (e.g. Bregman|1980). The most convincing kinematic 



information comes from Keeney et al. (2006) who found high velocity UV absorption along 
two sight-lines indicating a maximum outflow velocity of 250 km s _1 . 

We have discovered a population of compact, isolated clouds in the Australia Telescope 
Compact Array (ATCA) H I Galactic Center Survey that may trace neutral gas in a wind. 



2. A population of Clouds 



The ATCA H I Galactic Center Survey (McClure-Griffiths et al. 2012) covers Galactic 



longitudes —5° < I < +5°, and Galactic latitude —5° < b < +5° over the velocity range 
—309 km s _1 < Vlsr < +349 km s _1 with a velocity resolution of 1 km s _1 , an angular 
resolution of 145" and a mean rms H I brightness temperature of o"t 6 = 0.7 K. 

In the ATCA survey data for the central 1500 pc x 1500 pc of the Galaxy there is 
a population of compact H I clouds that are isolated both spatially and in velocity from 
the bulk of Galactic H I emission. Clouds were selected to be relatively compact (< 0.5°), 
discrete from other emission, and with brightness temperatures greater than 2<jT b over more 
than five spectral channels. Typically the clouds are only just spatially resolved. However, 
there are some more extended filamentary features, which we also refer to generally as 
clouds. Example clouds are shown in Figure [Tj We have found a total of 86 clouds between 
-209 km s- 1 < V LS r < 200 km s _1 . 

For each cloud we measure the central velocity, Vlsr, fitted peak brightness temperature 
after baseline removal, TJ,, velocity full-width half-max Av, mean radius, r, assuming a 
distance of R Q = 8.5 kpc, peak column density, Nnmax, and cloud mass, M c . The cloud 
properties are given in Table ll The mean radius, r = y/A/n, is calculated from the 
area of the cloud, A, where A is the number of pixels in the integrated intensity image 



with Njj > OANHmax- This is clearly a very rough estimate as many clouds are not well 
represented by an ellipse. The H I masses are calculated as the sum of all pixels greater 
than 40% Numax assuming a distance of 8.5 kpc. We quote H I masses assuming the gas 
is optically thin. The full width half max of the H I line (FWHM) is shown in Figure [2j 
along with histograms of mass, radius and column density for these clouds. The median, 
minimum and maximum values for the population are given in Table [2j There are no 
strong correlations of cloud radius, FWHM or mass with z height. The median line width 
corresponds to Tk = 22At> 2 < 4000 K. Assuming that the unresolved clouds are roughly 
spherical we estimate an average H I number density, (n c ) ~ Nu/(2r) = 1 cm -3 . 



There are several larger features, such as the loop at (7,6) = (359°, —3°) (Figure [Tj 
which appear to cohesive over many degrees. These features have also larger velocity widths 
and velocity gradients. The (/, b) = (359°, —3°) loop is particularly interesting because it 
appears to point directly back to the Galactic Center. 



2.1. Cloud Kinematics 

The location of the clouds, color-coded by their LSR velocities, is shown in Figure |3k,. 
The clouds clearly do not follow Galactic rotation. The velocity projection of an object at 
/, b onto the local standard of rest can be written 



V, 



LSR 



R sin I ( -£ - -^- ) - V R cos (I + 

ti JtCq 



cos b + V z sin b, 



where 9 is defined around the Galactic Center in a clockwise direction from the Sun-Galactic 
Center line. The velocity terms Vr, Vq and V z are the radial, azimuthal and vertical velocities. 
In circular Galactic rotation Vr and V z are zero and the sin I term produces a change in sign 
of Vlsr across I = 0°. The absence of this signature in the cloud distribution suggests that 
the kinematics of the clouds are not dominated by circular Galactic rotation. Because the 
clouds are at high velocities and do not show the signature of Galactic rotation we assume 
that the clouds are located at the Galactic Center. 

Clouds above the Galactic plane have predominantly positive velocities while clouds 
below the Galactic plane have predominantly negative velocities. Above the plane there are 
16 clouds at negative velocities and 29 clouds at positive velocities. Below the plane there 
are 28 clouds at negative velocities and 13 clouds at positive velocities. The mean LSR 
velocity at positive latitudes is 40.7 km s _1 and the mean velocity at negative latitudes is 
—52.8 km s _1 . This is unlikely to be the effect of the V z term in Equation 1, as sin b is very 
small and would require V z values of 5000 — 10000 km s _1 to explain the observed range of 
LSR velocities. 



This kinematic asymmetry is likely a product of confusion with unrelated H I emission, 
which effects our ability to find isolated clouds. It is well known that most of the H I emission 
around the Galactic center at vlsr > 100 km s _1 is at I > 0° and b < 0° and similarly the 
emission at vlsr < —100 km s _1 is at I < 0° and b > 0° (Burton & Liszt 1978 , see also 
Figure 4 of McClure-Grifnths et al. 2012). Consequently, we are not likely to find many 
high velocity clouds isolated in l-b-v space with positive velocities at b < 0° or with negative 
velocities at b > 0°. 



2.2. Simulated population 

To understand the selection effects on the distribution of observed cloud velocities we 
have simulated clouds in an outflowing wind. We have populated the interior of two cones 
having an opening angle, a, centered on the Galactic Center and extending to a height of 
1.5 kpc with 1000 randomly distributed clouds. Each cloud was assigned an LSR velocity 
based on its location, a constant conical wind of velocity, V w , and the expression similar to 



that derived by Keeney et al. (2006 



V, 



LSR 



v,„ 



(l — cos 2 b cos 2 /) =F cos b cos I 



P A 



(l — cos 2 b cos 2 I) 



1/2 



— 220 sin I cos b. 

(2) 

Here I and b are the Galactic coordinates, p is the distance along the line of sight to a cloud, 
r is the cloud's distance from the Galactic Center, and r 2 = Rq — 2pR & cos b cos/ + p 2 . The 
second term is negative for objects on the near side of the Galactic Center and positive for 
objects on the far side of the Galactic Center from the Sun. We assumed the IAU standard 
Rq = 8.5 kpc. 

Each simulated cloud was compared to the H I actually detected in the ATCA survey 
at its I and b. If the cloud lay at a velocity where there was no significant H I emission 
(< 4<ttJ it was assumed to have been "detected". Models covered wind velocities, V w , 
between 100 km s _1 and 800 km s _1 . Wind velocities < 150 km s _1 are unable to produce 
high velocity clouds as observed, whereas large wind velocities produce a cloud distribution 
whose mean velocity is higher than observed. Using a Kolmogorov-Smirnov (K-S) test, and 
separating the data into positive and negative latitude regions, we estimated the probability 
that the velocities of the observed clouds could be drawn from the same distribution as the 
simulated clouds. The K-S test p-values as a function of V w showed a peak near V w ~ 190 



lr This is a near reproduction of Equation 1 from 
the power of 1/2, correcting a typographical error in Keeney et al. 



Keeney et al. 



( 2006 1 except the second term is raised to 



km s _1 for positive latitudes and near V w ~ 220 km s _1 for negative latitudes, but at a lower 
confidence level. For all V w the simulated velocities were more consistent with the observed 
velocity distribution at positive latitudes than at negative latitudes. Adopting a threshold 
of p < 0.05, the K-S test implies that wind velocities greater than 270 km s _1 and less than 
150 km s _1 are not consistent with the observed cloud velocities. 

The observed cloud longitudes and latitudes set a lower limit on the cone opening angle, 
a. Values a < 3tt/4 radians are too narrow to produce intermediate latitude (\b\ ~ 2°) clouds 
at \l\ > 4°. The upper limit to a not well constrained as there is coupling between the angle 
and the assumed wind velocity. We experimented with several angles, including a fully 
spherical distribution, and found the results were best matched to a = 37r/4 radians. 

The resulting "detected" clouds for a wind velocity of V w = 190 km s _1 and a = 37r/4 are 
shown in Figure ^jp. Although the mean velocity of the initial distribution at both positive 
and negative latitudes was close to zero, because of the selection effects from blending with 
unrelated H I emission, the velocity distribution of the detected clouds shows the same 
preference for positive velocities at b > 0° and negative velocities at b < 0° observed in the 
real clouds. From 50 realizations of the simulation with V w = 190 km s" 1 we find mean 
and median positive latitude velocities of 52 km s _1 and 88 km s _1 , and mean and median 
negative latitude velocities of —4 km s _1 , and —60 km s -1 . The real data have mean and 
median values, respectively, of 40.7 km s _1 and 75.4 km s _1 at b > 0°, and —52.8 km s _1 
and —93.3 km s _1 at b < 0°. The velocity discrepancies at negative latitudes may indicate a 
higher wind velocity below the Galactic plane or a wind that does not point directly down. 

The simulation produces a large number of clouds at small longitudes (— 1° < I < +1°) 
and small velocities —70 km s _1 < Vlsr < +70 km s _1 at all latitudes. These clouds 
are collimated along the wind axis perpendicular to the Galactic plane. In most cases we 
expect these to be blended with unrelated low velocity H I emission. However, examining 
the observed H I data cube we do find significant structure collimated perpendicular to the 
Galactic plane at —1° < / < +1° with velocities |Vl5_r| ~ 50 — 70 km s _1 . These features do 
not meet our criteria of being discrete from other emission and were therefore not catalogued, 
but may reasonably be part of the same wind structure. 

The simulation thus approximately reproduces some of the key features of the observed 
distribution. The observations are consistent with a outflowing conical wind of ~ 200 km s _1 . 
We are also able to say with confidence that the observed velocity distribution of clouds with 
positive velocities at b > 0° and negative velocities at b < 0° is caused by confusion with 
unrelated H I. 



Discussion and Conclusions 



The recent discovery of the Fermi bubbles (Su et al. 2010) in 7-ray emission and their 



apparent counterparts in polarized radio continuum (Carretti et al. 2013) has fueled specu- 



lation that there is either a powerful jet from the Galactic Center or there has been a past 
burst of star formation. In explaining the polarized radio continuum structure of the Fermi 



bubbles, Carretti et al. (2013) argue that the bubbles are consistent with a past burst of 



star formation and confirm the Bland-Hawthorn & Cohen (2003) estimate that the bubbles 

10 55 erg. They estimate that the interiors of the Fermi bubbles 



required an injection of 

are filled with hot plasma of temperature T u 

traveling at a velocity of 1100 km s _1 . 



10 7 K with a density, n v 



9 x 10 3 cm 3 , 



The median radius, column density and velocity width of the Galactic Center H I are very 



similar to the lower halo tangent point clouds identified by Ford et al. (2010) and Lockman 



(2002). Ford et al. (2010) suggested that halo clouds are the remnants of supershells and 



superbubbles pushed beyond the Galactic disk through the stellar winds and supernovae of 
massive stars, but remaining dynamically linked to the Galactic disk. Our Galactic Center 



clouds have higher velocities than the Ford et al. (2010) clouds. The formation mechanism 



for the Galactic Center clouds, which determines their dynamics, is therefore probably more 
energetic than for the halo clouds. This is not surprising, the star formation rate in the 
inner 100 pc of the Galaxy is 0.1 M yr _1 (reviewed by Crocker 2012), whereas the star 
formation rate at the end of the Galactic bar where many halo clouds were found (I ~ 30°) 
is ~ 10~ 4 M yr" 1 (|Veneziani et al.||2013[). 



In spite of the extreme star formation rate in the very active region of the Galactic Center 
and estimated high temperatures in the wind, cold clouds are not unexpected. Simulations 



by Cooper et al. (2008) of starburst driven Galactic winds show compact clouds, filaments, 

10 cm~ 3 and temperature, T 



10 4 K entrained in 



and complex structures of density, n 
the inner 400 pc of a starburst wind. Similarly, the Melioli et al. (2013) simulations of the 
M82 wind produced a wealth of features at T ~ 10 4 K with sizes 20 — 300 pc and densities in 
the range 10 _1 — 10 cm -3 . The velocity of this cold phase of the outflow is typically 10 - 30% 



of the velocity of the hot phase, or 100 - 800 km s 1 , for the Cooper et al. (2008 ) and Melioli 



et al. (2013) models. Melioli et al. (2013) also found that in order for the wind to contain 



a rich filamentary structure of cool gas the starburst needed to occur over a relatively short 
timescale of a few million years. 

Our wind velocity modeling produces a velocity of the cold component of the wind that 



is similar to that observed in the Ha emitting filaments of M82 (Shopbell & Bland-Hawthorn 



1998) and predicted by hydrodynamic simulations (Cooper et al. 2008; Melioli et al. 2013). 



We can estimate the velocity of the hot component using the relationship for the ratio of the 



terminal velocity of the cold clouds, V t , to the hot wind, Vh w , for a cloud of density, n c , in a 



wind of density, n w , (Martin 2005) 



V t 



Vi 



liw 



3n u 
2n, 



11/2 



(3) 



Assuming V t = 

wind as Vhw ~ 



V w = 190 km s , as derived in §2.2 , we estimate the velocity of the hot 



1600 km s 1 , which is consistent with the hot wind velocity of 1100 km s 



estimated by Carretti et al. (2013) 



Presuming that the clouds are launched from very near the galactic plane, to reach 4° 
latitude in a ~ 200 km s _1 wind takes at least 3 Myr. If the clouds are entrained in a hot 
wind then they will be subject to a number of destructive effects. First, they will evaporate. 
The evaporation rate, rh, for an unmagnetized 15 pc cloud in a 10 7 K plasma with a density 
of n w ~ 9 x 10 -3 cm -3 is m ~ 1 x 10~ 3 M yr _1 , assuming saturated evaporation (Cowie 



et al.|l981 ). This implies a typical survival time of a few xlO 5 yr, which is less than the time 
for a cloud to reach a z height of 400 pc in a 200 km s _1 wind. However, the evaporation 
timescale is highly sensitive (oc T 5//2 ) to the assumed temperature of the plasma surrounding 
the clouds. Hydrodynamic models of winds show that the temperature may be anywhere 
in the range few x 10 6-7 K, it therefore seems possible that the clouds could survive for up 
to ~ 1 Myr. Also important to the clouds' lifetimes are the cloud crushing time, due to 
the wind shock, and the cloud cooling time. The relative timescales for these two effects 
gives an estimate of whether the cold, dense clouds would be protected by an enveloping 



radiative shock (Cooper et al. 2008). The cooling time for a cloud of density n c of shocked 
temperature T c ^ sh = 3/16 T w {n w /n c ) is on the order of a few hundred years. This is much 
less than the cloud crushing time, t cc ~ (n c /n w )(r c /v w ) ~ 4 Myr ( Cooper et al.|200"8 ), which 
suggests that a cloud could survive for some time in the hot wind. 

The total mass of H I in the observed clouds is only ~ 4 x 10 4 M . The mass loss rate, 
assuming the oldest cloud is ~ 4 Myr, is 0.01 M yr -1 . At this rate it is unlikely that this 
outflow will have a significant impact on the ~ 3 x 10 7 M reservoir of cold gas estimated for 



the 100 pc ring in the Central Molecular Zone (Molinari et al. 2011). We estimate the work 



required to accelerate these clouds to their current velocities; for a 270 M cloud traveling 
at 200 km s _1 , the kinetic energy is ~ 1 x 10 50 ergs, and the kinetic energy of the most 
massive cloud in our sample is ~ 8 x 10 50 ergs. The total estimated kinetic energy for all 
clouds traveling at a velocity of 200 km s _1 is ~ 1.5 x 10 52 ergs. Assuming the ratio of ratio 



of kinetic energy of the H I clouds to the total input wind energy is ~ 0.2 (Weaver et al. 



1977), as often assumed for superbubbles, the required wind energy input was 

ergs over the past 10 6 yrs. The wind luminosity is therefore ~ 2 x 10 39 erg s" 

easily supplied through supernovae and stellar winds in the GC region with a star formation 



~ 8 x 10 52 
l , which is 



rate of ~ 0.1 M Q yr" 1 (|Crocker|[20T2|. 



4. Summary 

We have identified a population of compact, cool, T < 4000 K, H I clouds around 
the Galactic Center. Although physically similar in terms of size, temperature and column 
density to other cold H I clouds in the lower halo of the Galaxy, the Galactic Center H I 
clouds have unique kinematics and do not follow Galactic rotation. We suggest that these 
clouds are associated with a wind from the Galactic Center and are the remnants of cold, 
dense gas entrained in the outflowing wind. Using a simple model for a Galactic wind with a 
velocity of ~ 200 km s _1 we are able to roughly reproduce the observed velocity distribution 
of the clouds. This value, and hydrodynamic simulations of a starburst wind, which find 



that the velocity of the cold phase gas is 10 - 30% of the velocity of the hot phase ( Cooper 



et al. 2008; Melioli et al. 2013), the hot wind emanating from the Milky Way's Center must 
have a velocity of ~ 600 — 1900 km s _1 . We estimate the luminosity required to drive the 
clouds is ~ 2 x 10 39 erg s _1 , which is easily supplied by stellar winds and supernovae in the 
Galactic Center region over the past ~ 2 x 10 6 years. 

We are grateful to R. Crocker for useful comments on an early draft of this manuscript. 
The ATCA is part of the Australia Telescope which is funded by the Commonwealth of 
Australia for operation as a National Facility managed by CSIRO. 
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Fig. 1. — H I images at (a) v = 161.5 km s 1 , (b) v = 83.2 km s 1 , (c) 



■121.2 km s~\ 



and v = —90.7 km s . Catalogued clouds whose central velocity lies within 5 km s of the 
shown velocity are circled. 
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Fig. 2. — Histograms of: (a) the mean radius of individual clouds, (b) the velocity FWHM 
of the clouds, (c) the peak column density, Nh, of the clouds and (d) the total mass of the 
clouds. 
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Fig. 3. — Distribution of clouds around the Galactic Center found in (a): the ATCA Galactic 
Center survey and (b): a simulated wind of velocity 190 km s . The color of the symbols is 
determined by the central velocity of the sources. The LSR velocity of the simulated clouds 
is calculated from Equation [2j 
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Table 1. Compact clouds catalogue 



Name 



I 

(deg) 



b 
(deg) 



Vlsr 

(km s -1 ) 



Peak T b 

(K) 



Av 
(km s _1 ) 



Nh max 

'xlO 19 cm" 



M c 
(M ) 



(PC) 



GO. 2+3.8+65 


0.23 


3.81 


63.7 


4.4 


5.9 


5.7 


25 


6.6 


GO. 2-1.6-170 


0.27 


-1.62 


-169.6 


4.9 


17.5 


17.6 


2190 


32.9 


GO. 3+2. 8+195 


0.30 


2.84 


200.8 


1.4 


28.6 


9.7 


220 


12.2 


GO. 3+3.7-82 


0.36 


3.74 


-74.3 


3.7 


13.7 


11.1 


706 


22.4 


GO. 7+4.5+61 


0.68 


4.52 


63.3 


9.6 


5.3 


10.8 


287 


14.2 



Note. — Table [T] is published in its entirety in the electronic edition of the Astrophysical Journal Letters. 
A portion is shown here for guidance regarding its form and content. 



Table 2. Cloud population properties 



Property 



Median Min 



Max 



r (pc) 15.0 


3.5 


32.9 


Av (km s" 1 ) 13.6 


3.0 


30.9 


T b (K) 3.8 


1.4 


14.2 


N Hmax (10 19 cm" 2 ) 9.9 


1.0 


34.9 


M c (M ) 270 


5 


2.1 x 10 3 



